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Abstract. We present a solution to obtain a high-resolution image of a wide field with the central source removed 
by destructive interference. The wide-field image is created by aperture synthesis with a rotating sparse array 
of telescopes in space. Nulling of the central source is achieved using a phase-mask coronagraph. The full (u,v) 
plane coverage delivered by the 60m, six 3-meter telescope array is particularly well-suited for the detection and 
characterization of exoplanets in the infrared (DARWIN and Terrestrial Planet Finder (TPF) missions) as well 
as for other generic science observations. Detection (S/N=10) of an Earth-like planet is achieved in less than 10 
hours with a Ifim bandwidth at 10/^m. 

Key words. Techniques; interferometric - Instrumentation: high angular resolution 



1. Introduction 

As seen from 10 parsec ( pc) away, the Earth is 0.1 arc- 
second (") away from the Sun and 10^ times fainter at 
visible wavelengths. At 10 fj,m, the light intensity ratio is 
only 10^, and a 20m minimum baseline is needed to resolve 
the system. The lower Planet /Star intensity ratio as well 
as the likely presence of biomarkers (CO2, O3, i?20, CH4) 
in the atmospheric spectrum from 5 /zm to 15 /zm led us to 
choose this wavelength range for our study. A space inter- 
ferometer with nulling capabilities is well suited for this 
work because of the large baseline needed to resolve the 
system and the high luminosity ratio between the planet 
and the star. 

Bracewell (1978) proposed to combine the beams 
from 2 telescopes to form interferometric fringes with 
the star centered on a dark fringe (optical path- 
length difference of half a wavelength between the 2 
light beams). The on-axis star's light would then be 
nulled, and by rotating the 2-telescope interferome- 
ter around its optical axis, the light from a compan- 
ion would be temporally modulated. The same con- 
cept has more recently been generalized to several tele- 
scopes to obtain a deeper null and a better (u,v) cover- 
age dAngel fc Woolf 1994 |Menncsson fc Mariotti 1991\j. 



A monopixel detector records a light intensity, which is 
the integral of the product of the light distribution on 
the sky and a "transmission map" given by the geometry 
of the array. As the interferometer rotates, the residual 
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light of the on-axis star remains constant and an off-axis 
companion is revealed as a temporal modulation of the 
recorded signal. The imaging capabilities of these concepts 
are very limited (small field of view and very partial (u,v) 
plane coverage) , increasing the risk of confusion between 
an exoplanet and an anisotropy of the exozodiacal cloud. 
The detection by light modulation on a monopixel detec- 
tor suffers from high sensitivity to temporal variations of 
both the thermal background intensity and the light leak- 
age from the nulled central star. In this paper, we present 
a solution to image (with a full (u,v) plane coverage up 
to the cutoff frequency) a field with the central source re- 
moved by destructive interference. The wide field of view 
and good (u,v) plane coverage limit the risk of confusion 
with background sources and exozodiacal anisotropies. We 
demonstrate the importance of good (u,v) coverage and 
show how to achieve it in §2. The pupil densification and 
pupil redilution techniques are presented in §3: these tech- 
niques make possible the use of a nulling coronagraph on 
a sparse array of telescopes. The overall behavior of the 
concept is studied in §4. In §5 and §6, we explore the ca- 
pabilities of this concept to image Earth-like planets and 
study other astrophysical objects. 



2. (u,v) plane coverage 

2.1. Importance of good (u,v) plane coverage 

Very partial (u,v) plane coverage is sufficient for observa- 
tion of objects for which a simple model exists. Ground- 
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based interferometers are routinely measuring stellar di- 
ameters using 2 telescope-interferometers. By changing 
the projected baseline (by physically moving the tele- 
scopes or taking advantage of Earth's rotation), more 
points are accessible in the (u,v) plane and limb darkening 
can be measured (Bajian et al. 1998). Unfortunately, no 
simple model exists for a planetary system, and exoplanet 
imaging missions like TPF (Beichman 1998) and Darwin 
(Leger & Mariotti 1996) need good (u,v) plane sampling 
before detecting with certainty an exoplanet. A planetary 
system is likely to have several planets embedded in a zo- 
diacal cloud. As seen face-on, a planetary system with an 
Earth- like planet 1 Astronomical Unit {AU) from a Sun- 
like star in a 1 Zodi (unit of dust density in a zodiacal 
cloud, normalized to the solar system zodiacal cloud) zo- 
diacal cloud is as bright as a 0.3 AU by 0.3 AU "pixel" of 
exozodiacal light next to it. 

The presence of massive bodies in the zodiacal cloud 
will perturb the orbits of dust particles and produce an 
elongated concentration of zodiacal dust: in the solar 
system, the Earth is responsible for a 10% variation of 
the surface brightness of the zodiacal dust as seen from 
Earth (Reach ct al. 1995, Dcrmott ct al. 1994). Because 
in some cases these clumps can be brighter than the 
planet itself ( Liou fc Zook 199^ ), it has been proposed 
to infer the existence of a planet from observation of 
the exozodiacal cloud. Recent direct ground-based ob- 
servation of the dust stream responsible for the Leonid 



meteor shower (Nakamura et al. 200C) and IRAS/ISO 
observations of cometary trails (pykes fc Walker 199^ , 



Kresak 1993, Jenniskens et al. 1997) demonstrate the 



long-lasting effect of cometary debris. In our solar sys- 
tem (seen face-on), the light intensity from dust concen- 
trations in the zodiacal emission (excluding the smooth 
component of the zodiacal cloud) , when integrated over a 
scale of a tenth of an AU, is less than a percent of the 
Earth's light intensity. However, in systems with a higher 
zodiacal component or no massive planet outside of the 
Earth's orbit, dust concentrations could be as luminous 
as an Earth-like planet. This becomes even more critical 
if the system is seen edge-on, since the surface-brightness 
of planet-induced arcs of dust would be amplified when 
seen tangentially. 

The TPF and Darwin missions propose to image stars 
at a distance up to 20 pc to look for planets. Being able to 
rule out arcs of exozodiacal light excess and background 
sources will make the detection of planets more reliable 
and, to some extent, reduce the need for repetitive ob- 
servations of the same object to confirm a detection. A 
large field of view is also essential for characterization 
of solar systems from 5 to 20 pc. While imaging a Sun- 
Earth system at 20pc requires at least a 50 milli-arc-second 
( mas) resolution (separation corresponds to 1 resolution 
elements or more) , the angular separation between a Sun- 
Jupiter system at 5 pc is 1 " . A wide field of view requires 
a well-filled (u,v) plane, which also allows observations of 
various astrophysical sources for which no reliable model 
exists. 



2.2. Obtaining full (u,v) plane coverage with aperture 
synthesis 

Several technical constraints limit the achievable (u,v) 
plane coverage: number of telescopes, telescope diame- 
ter and array geometry. A large number of telescopes or 
a larger diameter of individual telescopes would increase 
the cost and complexity of the mission. As the number of 
telescopes increases and their diameter decreases, the in- 
ternal metrology of the array becomes more complex and 
less photons per aperture are available, making precise 
high-speed fringe tracking more difficult. We have chosen 
to work at 10 /im, and to reach the resolution needed to 
separate a Sun-Earth system at 20 pc at this wavelength, 
a 60 TO baseline is needed (34TOas resolution at 10/ito). To 
further limit the number and size of telescopes, we have 
chosen to use rotation of the whole array to fill the (u,v) 
plane. Rotation of the array is the easiest way to change 
the set of (u,v) points measured, since it does not require a 
change in distances or relative angles between apertures (a 
rigid array can be used) and there is no need for large opti- 
cal pathlength differences corrections. Rotation of a rigid 
interferometric array greatly improves its (u,v) plane cov- 
erage. By adopting an optimal configuration, it is possible 
to obtain a full (u,v) plane coverage through rotation of 
an array of small telescopes. Guyon and Roddier (2001) 
have explored this problem for arrays of up to 10 identical 
telescopes and, for a given baseline, computed the con- 
figurations that allow full (u,v) plane coverages with the 
smallest possible telescope diameters. Table 1 lists, for a 
given number of telescopes, the minimal telescope diam- 
eter required to fill the (u,v) plane up to the resolution 
limit of the array (60 m) . 



N 


4 


5 


6 


7 


8 


9 


10 


D(m) 


4.62 


2.90 


1.98 


1.43 


1.11 


0.91 


0.76 



Table 1. Minimum telescope diameter (D) to fully 
cover the (u,v) plane up to the maximum frequency 
for a 60m array of N telescopes. Adapted from 



(Guyon fc Roddier 2001) 



The array configuration for each of those options is 
given by Guyon & Roddier (2001). From this table, six 
2 TO telescopes can fully cover the (u,v) plane up to the 
maximum frequency of a 60m array. In this work, when- 
ever numerical values are given, we choose to adopt the 
6-telescope configuration, with d = 3 to (diameter of each 
aperture of the interferometer) to obtain a good SNR at 
all spatial frequencies. The array geometry is shown in fig. 
1 (upper left). In this paper, N is the number of apertures, 
d the diameter of each aperture and A is the wavelength. 
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Fig. 1. Array configurations (top) and corresponding 
PSFs (bottom) for a non-redundant (left) and a redun- 
dant (right) array geometry. The diffraction peaks in the 
PSF of the redundant array are more contrasted. The non- 
redundant six telescope array on the top left is optimized 
for full (u,v) plane coverage when combined with rota- 
tional aperture synthesis. 



2.3. Tradeoff between (u,v) plane coverage and planet 
detection sensitivity 

Obtaining a complete (u,v) plane coverage with a limited 
number of small telescopes requires a non-periodic array. 
The Point Spread Function (PSF) obtained from a non- 
redundant array is poorly contrasted whereas a highly re- 
dundant array yields very contrasted series of diffraction 
peaks (Fig. 1). For planet detection, our concept uses both 
nulling (canceling most of the starlight photons) and imag- 
ing (separating the photons of the planet from the pho- 
tons of the star on the detector). The lower contrast of the 
diffraction peaks reduces the signal to noise ratio (S/N) 
for detection of exoplanets because the spatial separation 
of the planet light and the star light is less efficient. We 
give some estimate of this decrease in S/N in Appendix 
B. 

The efficiency of the (u,v) plane coverage could be sac- 
rificed to gain S/N for planet detection at the expense of 
a reduced field of view and lower imaging capabilities (in- 
creased chances of confusion). The poorer imaging capa- 
bilities translate into a lower confidence in the detection, 
requiring multiple observations (longer exposure time). It 
is hard to estimate the overhead of a poorer (u,v) plane 
coverage because it is very dependent upon the structure 
of the object observed. For observations of complex astro- 
physical sources (galaxies, star forming regions) for which 
nulling is not critical, there; is no such tradeoff and the 
coverage of the (u,v) plane should be optimized regardless 
of the instantaneous PSF contrast. 



The nulling performance of our concept is independent 
of the array geometry. It is therefore possible to adopt a 
redundant configuration without any change in the con- 
cept. A reconfigurable array is also a possible solution for 
choosing the best array geometry for each type of obser- 
vation. 

In this work, we have adopted a nonredundant con- 
figuration which obtimizes the (u,v) plane coverage. This 
makes our conccipt a true imager, with excellent scien- 
tific capabilities for projects others than exoplanet detec- 
tions. The risk of confusion between exoplanets, exozodi- 
acal structures and background sources is also very low. 
Alternatively, the redundant array shown in Fig. 1 (top 
right) might also be considered for exoplanet detection 
and characterization: the better contrasted PSF allows 
this array to reach the same S/N ratio with exposure times 
3 to 4 times shorter (see Appendix B). Although such a 
choice would be risky for observation of complex targets, 
this array would be most efficient on systems with a small 
number of planets and little exozodiacal structure. 



2.4. The image reconstruction algorithm 

A scries of "snapshot" images is acquired with the array, 
with a direct recombinaition of the beams at the focal 
plane (Fizeau imaging). Each snapshot is acquired at a 
given rotation angle of the array. An example of aperture 
synthesis using this array, with rotation, is given in Guyon 
& Roddier (2001). An important step of the reconstruction 
is the Fourier filtering of each snapshot frame to reduce the 
photon noise : Fourier components that are not sampled 
by the current configuration of the array are set to zero 
for each frame. 

The rotation of the array can be continuous, provided 

that, if the snapshot exposure times are long (the array 
rotates significantly during one snapshot), the focal plane 
detector is counter-rotating in order to maintain a fixed 
orientation relative to the sky. Thanks to the linearity 
of the Fourier transform, each snapshot then samples a 
set of Fourier-space domains which are thick arcs rather 
than spots : the position of the objects on the detector is 
constant during an exposure, but the changing (rotating) 
diffraction patterns are integrated through the exposure. 
Each snapshot is equivalent to an image acquired by a 
60m telescope with mask in the pupil plane having six 3 m 
diameter holes. Fizeau imaging is preferred to the classical 
fringe amplitude and phase measurements (one fringe per 
pair of telescopes) because collapsing the beam of each 
telescope into a coherent wave would limit the total field 
of view to X/d (0.7 " from edge to edge) which, in this 
example, is too small to image a Sun-Jupiter system at 
10 pc. 



3. Pupil densification and nulling coronagraphy 
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Fig. 2. Cut through a diluted pupil's PSF (top) and a 
densified pupil's PSF (bottom). A 2 apertures interferom- 
eter was simulated to generate the 2 plots. The envelope 
(PSF of a single aperture) is over-plotted on the PSF in 
each case. Because the distance between the 2 apertures 
is smaller when the pupil is densified, the fringe period is 
larger. 




Fig. 3. Pupil (left) and PSF (right) in the densified pupil 
configuration for a 6 telescope array. 



3.1. Introduction to pupil densification 

The PSF of an interferometer (of baseline B and sub- 
apertures diameter d) in Fizeau imaging is characterized 
by two quantities : 

— The size of the central diffraction peak (A/i?) 

— The size of the PSF's envelope in which the diffraction 
peaks he (A/d). 



Pupil densification increases the pupil filling factor by 
bringing the sub-pupils closer to each other, therefore 
artificially increasing d/B and reducing the number of 
diffraction peaks in the PSF. It "magnifies" the diffrac- 
tion pattern inside the envelope (fig. 2). When the pupil 
is fully densified, only one bright diffraction peak is inside 
the PSF's envelope. While the PSF is field-invariant in 
the diluted pupil scheme, it is not in the densified pupil 
scheme. Pupil densification of a sparse array of apertures 
creates a PSF close to a single-aperture telescope's PSF 
(Labeyrie 1996, Pedretti et al. 2000) which allows a phase 
mask coronagraph to work efficiently. Figure 3 shows the 
densified pupil and corresponding on-axis PSF and is to be 
compared with the PSF in the diluted pupil configuration 
(Fig. 1, bottom left). In the diluted pupil configuration, 
the PSF has many bright secondary peaks while in the 
densified pupil configuration, most of the light is concen- 
trated in the central peak. 



3.2. Requirements for coronagraphy and wide field of 
view imaging 

Because of the non-redundancy of the array (good (u,v) 
coverage), the diffraction peaks of the PSF in the diluted 
pupil scheme tend to break into multiple fringes, with 
lower contrast, away from the center of the PSF (Fig. 1). 
If a coronagraph were to be used in this focal plane, it 
would have to "mask" each of those peaks. Because of the 
"fading" of those peaks into multiple low-contrast fringes, 
this "masking" would also mask most of the light of any 
point source in the field of view. Another problem is the 
wavelength-dependence of the PSF. Diffraction peaks of 
the PSF are spectrally dispersed in a direction pointing 
to the center of the PSF. The combination of these two 
effects makes it impossible to achieve a good nulling of the 
star without also masking the planet. 

A solution to the first problem would be to make 
the array redundant (regular spacing between aper- 
tures) to form high-contrast diffraction peaks across the 
PSF. The wavelength effect can be solved by using a 
wavelength-dependent magnification device to cancel out 



the wavelength-dependent scale of the PSF (Wyne 1979 



Roddier et al. 1980). Another solution to the wavelength 



problem is to use a ID array to free up one dimension 
in the focal plane: this dimension can then be used to 
spectrally disperse the PSF on a coronagraphic device 
which follows the changing scale of the PSF on this axis. 
This is the solution that Aime et al. (2001) have explored 
with a linear periodic array. With a non-redundant, (u,v) 
coverage-optimized, array, pupil densification is needed to 
concentrate the light of the central source into a single 
high-contrast diffraction peak on which coronagraphy can 
be performed: the entrance pupil of the array is densified 
into a tight configuration before coronagraphy is applied. 

Formation of a wide field of view {FOV) image 
(FOV > X/d) requires Fizeau- mode imaging with a non- 
densified exit pupil. The number of resolution elements in 
the "clean field of view" of the densified pupil is too small 
to recover the information needed to create a wide field 
of view image using rotational aperture synthesis. This 
is why, after the coronagraphic device, the pupil needs to 
be rediluted into its original configuration (entrance pupil) 
before Fizeau combinaition of the beams on the focal plane 
detector (Fig. 4). Table 2 gives the Field of view and the 
number of diffraction peaks in each pupil configuration. 

3.3. Sub-Aperture configuration in the densified pupil 

Coronagraphy works best when the densified pupil is tight 
(high filling factor). Such a tight configuration cannot be 
obtained if the relative positions of the sub-pupils' centers 
are kept constant: the pupil needs to be rearranged for an 
efficient densification. This seems to violate a common law 
of image formation through interferometers which states 
that the relative position of the sub-pupils of an interfer- 
ometer should not be modified before image formation. In 
fact, this rule is not violated because the final image will 
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Densified 


Diluted 




Pupil 


Pupil 


Field of 


Useful FOV 


No 


View 




limitation 


Number of 


1 


^ ds/N' 


diffraction peaks 



Table 2. Field of view and number of diffraction peaks of 
the PSF in the densified pupil and diluted pupil schemes. 
N is the number of apertures, B the baseline and d the 
diameter of individual apertures. 




liitci"fcromctcr Imaging fi 
exit pupil plane 



Fig. 4. The entrance pupil of the interferometer is den- 
sified before the coronagraph. The densified exit pupil of 
the coronagraph is rediluted before imaging on the detec- 
tor array in order to maintain a wide field of view. 

be created after redilution of the array into its original 
configuration. 

However, this means that the light complex amplitude 
at the coronagraphic focal plane does not follow the sim- 
ple laws of the formation of an image. For an on-axis 
point source, the wavefront phase is constant on all sub- 
apertures and across them. This flat wavefront in the den- 
sified pupil will yield a sharp PSF in the coronagraphic 
focal plane. If the point source is not on-axis, wavefront 
slopes are present on each sub-aperture and large phase 
shifts exist between sub-apertures. In the entrance pupil, 
those wavefronts are all part of a large tilted wavefront, 
and the tilt corresponds to a translation of the PSF. In 
the densified pupil, the slopes and phase shifts cannot be 
obtained by multiplication of a tilted wavefront by the 
densified pupil function (1 inside the pupil and outside): 
the PSF is not translation-invariant. For some positions 
of the point source, the PSF in the densified pupil case 
will be broken into several low-contrast diffraction peaks, 
while for others, the set of fringes will coherently add up 
into a single diffraction peak (Fig. 5). For an on-axis point 
source, the coronagraph will "mask" the light of the source 
because a clear diffraction peak is centered on the optical 
axis. 
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Fig. 5. PSFs in the coronagraph's focal plane (after pupil 
densification) for different positions on the sky. The lower 
left PSF corresponds to an on-axis point source. In this 
grid, the point source position in the sky is incremen- 
tally increased by lOmas steps along two orthogonal axis. 
The white circle plotted on the on-axis PSFs indicates the 
"clean field of view" in the densified pupil scheme. 



3.4. Nulling coronagraphy 

Pupil densification reduces the useful field of view in the 
coronagraph's focal plane to about ViV x times a/ZV x . 
With a small number of telescopes {N < 10), this field 
of view is only a few times bigger than the diffraction 
peak of the PSF. There is therefore a very strong con- 
straint on the size of the mask to use in this concept: 
a coronagraph whose mask would "block" an area equal 
to or larger than \/N x ^ times Vn x ^ in the coro- 
nagraph's focal plane would in fact "block" the light of 
any source closer than ^ to the optical axis. Figure 6 il- 
lustrates this effect by showing the transmission map of 
the concept on the sky when a phase mask coronagraph is 
used: the distance between consecutive nulls is comparable 
to the width of each null. If < 10, a Lyot coronagraph 
shouldn't be used since a Lyot coronagraph, even with an 
apodized mask, cannot reach good null performance (bet- 
ter that 10'^) with masks smalle r than 2 x The phase 



mask (Roddier & Roddier 1997, Guyon et al. 199!;), how- 



ever, offers a very attractive solution because of its very 
small size (less than half the diameter of the first dark 
Airy ring for a circular aperture) and its theoretical to- 
tal extinction for an on-axis point source. The four quad- 



rants phase mask (Rouan et al. 2000, Riaud et al. 2001) 



offers similar nulling performance but would yield a lower 
overall transmission due to partial extinction on the two 
perpendicular axes of the phase shifts, and is sensitive 
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to the shape of the densified pupil (the four quadrants 
phase mask would not work on the densified pupil ge- 
ometry wc have adopted in this work). The phase mask 
coronagraph can work efHciently with almost any pupil 
shape, provided that the pupil is not too sparse, but the 
apodization mask required in the entrance pupil of the 
coronagraph (Appendix A) is different for different pupil 
shapes: if a single apodization mask is used, the densified 
pupil shape cannot be modified. However, the four quad- 
rants phase mask does not suffer from the wavelength- 
dependant PSF scale. For N < 10, the lower overall trans- 
mission of the the four quadrants coronagraph seriously 
affects the sensitivity of the study. However, for large val- 
ues of N, this becomes a less serious problem, and the 
advantages of the four quadrants phase mask (no need 
for a pupil apodization mask, insensitivity to wavelength- 
dependant PSF scale) could make it a more attractive 
solution than the phase mask. 

Therefore, we have chosen to use the phase mask coro- 
nagraph with an apodized densified pupil. This technique 
offers a total extinction of an on-axis monochromatic point 
source. Details of the apodization technique as well as per- 
formance and sensitivity of the phase mask to various er- 
rors are given in appendix A. The results of this detailed 
study of the phase mask coronagraph technique are used 
to estimate the sensitivity of the concept to point-source 
detection (Appendix B). 

4. Imaging with tlie pliase mask coronagrapli 
applied on a sparse array of telescopes 

4.1. PSF characteristics 

Introducing a coronagraph after densification of the pupil 
(and before rcdilution) removes the invariance by transla- 
tion of the PSF in the final focal plane (after pupil redi- 
lution). The light of an on-axis point source is stopped by 
the coronagraphic assembly while for other point source 
positions, most of the incoming light is detected in the 
final focal plane. This field-dependent modulation of the 
transmission (Fig. 6) is accompanied by a change of the 
PSF structure (Fig. 7). 

4.2. Array rotation and residual star light removal 

The image of a star (disk of uniform brightness) on the 
focal plane is invariant with the rotation of the array (the 
detector is corotating with the array). The image of an 
off-axis point source (planet), however, is rotating around 
the optical axis and changing in structure (effects intro- 
duced by the coronagraph, see §4.1). It is hence possible 
to decouple the static image of the star (and other rota- 
tionally symmetric light distributions) from the changing 
image of planets (and other non-rotationally symmetric 
light distributions). For example, the median of all the 
frames gives a very good approximation of the residual 
star light and can then be subtracted from each indi- 
vidual frame. When a high number of snapshots have 



V 




Fig. 6. Transmission map of the interferometer on the sky. 
The scale is linear from (black, no transmission of the 
source's light) to 1 (white, total transmission). The image 
is 1.024" by 1.024". 




Fig. 7. PSFs of the concept for 4 different positions on 
the sky. Each frame is 1.6" by 1.6". Positions on the sky, 

in mas, relative to the optical axis, arc (+10, +10) (upper 
left), (+20,+20) (upper right), (+30,+30) (lower left) and 
(+40,+40) (lower right). 



been acquired at various position angles of the interfer- 
ometer, the median of all snapshots yields a very reli- 
able image of the rotationally symmetric component of 
the observed object (Star + rotationally symmetric com- 
ponent of the exozodiacal cloud). For detection of exoplan- 
ets (or other non-rotationally symmetric structures), this 
self-calibrating technique is expected to be significantly 
more reliable than observation of reference stars. Thanks 
to this technique, any noise which is not correlated with 
the position angle of the array (such as "speckle" noise) 
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will be efficiently decoupled from the planet's signal and 
will average out efficiently. 

Our simulations show that the distribution of the resid- 
ual starlight on the detector has little dependence with 
the angular size of the star (provided that the star is not 
resolved) . It can be approximated by a reference distribu- 
tion multiplied by a coefficient which increases with the 
angular diameter and surface brightness of the star. This 
property allows precise modeling of the residual starlight 
in the focal plane. Thanks to the very large number of 
measurements available in an image (as opposed to a 
monopixel detector) each frame contains significant in- 
formation about the status of the interferometer. Slight 
pointing errors of the array (equivalent to a phase error 
between the sub-pupils) can be diagnosed through analy- 
sis of the residual starlight distribution. The frames can 
be used in a low frequency feedback loop to correct for the 
non-common path phase delays between the fringe track- 
ers and the detector array: the fringe tracker would intro- 
duce a constant phase delay between the beams to correct 
for these delays. 

To remove the residual star light for detection of exo- 
planets, subtraction of the rotation-invariant component 
of the light distribution on the detector is a robust ap- 
proach. Thanks to the high information content of each 
frame, small variations of the central star image can 
be modelled in order to further improve this technique. 
Stellar disk elongation (for rapidly rotating stars seen 
"edge-on") and variations of the null are the two major 
potential causes of such variations. 

4.3. Final image reconstruction 

A final image can be created from the snapshots using 
rotational aperture synthesis ( puyon fc Roddier 2001 ). 
However, the presence of a coronagraphic device degrades 
the quality of the reconstruction: the Fourier transform of 
a frame does does not yield "pure" Fourier components of 
the real object image. The frame is affected by the coro- 
nagraph and no simple exact relation exists between the 
frame and the object. This "degradation" of the recon- 
structed image becomes more serious as the number of 
telescopes is reduced. We believe it becomes acceptable 
for an array of 6 telescopes or more: for a 6 telescopes ar- 
ray, this effect decreases the S/N by a factor of 3.5 (§5.3). 
Although this factor is quite large, for arrays of 6 or more 
telescopes, the efficient use of the photons and good imag- 
ing capabilities makes this concept very attractive when 
compared to "monopixel nuUers" . Since the formation of 
the image is a linear process with a field-dependent PSF 
it should be possible to correct for this effect. However, 
the construction of a deconvolution algorithm specific to 
this problem has not been explored in this work. Figure 8 
shows a reconstructed image illustrating this effect. 



Fig. 8. Simulated image of the solar system's planets at 
a distance of 10 pc. The values adopted for angular sepa- 
rations and fluxes are given in Table 3. No photon noise 
has been included. The image of the Sun was numerically 
substracted. 

5.1. Image reconstruction for complex planetary 
systems 

In this section, we illustrate the potential of our concept at 
recovering images of complex planetary systems by using 
a non-redundant array of 6 apertures optimized for (u,v) 
plane coverage. Rotation of the array is used to fully cover 
the (u,v) plane. We simulated the observation of a solar 
system at a distance of lOpc. Zodiacal light and exozodia- 
cal light have not been included in this simulation. Photon 
noise was not taken into account (infinite exposure time). 
The characteristic values for the planets in this model are 
given in Table 3. 





6(mas) 


Teff(K) 


Fjy 


Fph 


Sun 







3.2 


5 10"" 


Mercury 


39 


480 


5.8 10"' 


0.91 


Venus 


72 


250 


2.2 10"' 


0.35 


Earth 


100 


300 


6.4 10"' 


1.0 


Mars 


152 


270 


1.0 10"' 


0.16 


Jupiter 


520 


140 


3.3 10"' 


0.52 


Saturn 


950 


110 


1.4 10"** 


0.022 


Uranus 


1914 


68 


7.7 10"'^ 


1.2 10"" 


Neptune 


3000 


55 


5.0 10""" 


7.7 lO""* 



Table 3. Values adopted in the simulation for the angular 
separation 9 to the sun, effective surface temperature (at 
lOfim) and flux {Fjy in Jy and Fph in p/i.s"^.m"^./im"^) 
for a solar system as seen from 10 pc away. 



5. Planet detection capabilities 



Figure 9 shows the (u,v) plane coverage of the obser- 
vation (half a rotation of the array). The reconstructed 
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Fig. 9. Map of the (u,v) coverage obtained after half a 
rotation of the array. A cut along the diameter of this 
disk is also shown (linear scale). 

image, in which the image of the sun has been subtracted, 
(Fig. 8) demonstrates the ability of this concept to im- 
age planets close to the null (Mercury) as well as planets 
further out (Saturn). The good quality of the image re- 
construction over a wide range of distances allows unam- 
biguous detection of the 6 brightest planets (Uranus and 
Neptune are very faint due to low effective temperature 
at 10/im). Arcs of exozodiacal light and other resolved 
structures could also be imaged accurately. The residual 
"cloud" of speckles is due to the the effect of the corona- 
graph, which corrupts the (u,v) plane values obtained on 
each snapshot. We have not attempted to remove those 
predictable artifacts from the image by deconvolution. 
Because we simulated a system with no exozodiacal light, 
the effect of this non-perfect reconstruction algorithm is 
easily visible, since, without the use of a coronagraph (and 
without a central star), the background between the plan- 
ets should be null. 

5.2. Theoritical detection sensitivity in imaging mode 

In Appendix B, we give an theoritical estimate of the 
point source detection sensitivity in a snapshot frame. In 
this paragraph, we use this estimate to quantify the point 
source detection sensitivity in imaging mode. The estimate 
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Transmission in imaging mode 




Instance from optical axis (mas) 

Fig. 10. Transmission of the interferometer as a function 
of distance to the optical axis in the imaging mode. 

in Appendix B was computed for a rotation angle of the 
interferometer such that the planet's image is not occulted 
by the coronagraph. In a real imaging observation (series 
of snapshots as the interferometer rotates), the planet's 
image is periodically occulted by the coronagraph. While 
on some snapshots the S/N of the detection is given by 
equation B.23, on others, the planet's signal is partially 
rejected by the coronagraph and the S/N of the detection 
is lower. Figure 10 shows the average transmission, T, of 
the interferometer as a function of distance to the optical 
axis, X. At a distance of 0.1", the average transmission is 
50%. On the final reconstructed image, the planet's flux 
is the sum of the planet's flux in each snapshot. Following 
the notations adopted in Appendix B, since this average 
transmission factor affects Npi, Nz and Nez, but not Ngti 
and Nst2, the expression for the signal to noise of a point 
source detection becomes 



N ^ ^/{Npi + Nz + Nez)T{x) + N,ti + Nst2 

With T(.x) = 0.5 and the values computed in Appendix 
B, the exposure time required to reach a given S/N is 3.67 
times longer than the exposure time required to reach the 
same S/N in a snapshot in which the planet is not occulted 
by the coronagraph. Therefore, in imaging mode, with the 
non-redundant six 3 m apertures, detection {S/N = 3) of 
the Earth at IQpc is achieved in 2Qmn.ixm {hh.jim with an 
0.1 efficiency). In this work, exposure times will be given 
in time x bandwidth units, such as s./xm, corresponding 
to a 1/xm bandwidth. Is./xm corresponds to a Is exposure 
time for a 1/im bandwidth, or a 0.2s exposure time for a 
2^m bandwidth. 

T{x) is higher, except in the central null, when the 
number of telescopes is larger. 

5.3. Numerical simulation with photon noise 

For a solar system at lOpc, we have identified in Appendix 
B that the major source of noise is the photon noise from 
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the main star residual light (due to partial resolving of its 
disk). The effect of "speckle" noise, arising from variations 

in the wavcfront of individual beams and OPD variations 
(this definition includes pointing errors) will be kept small 

— With a wavefront accuracy of 25nm, the dominant 
source of light contamination arises from the partial 
resolving of the star's disk (Appendix B). Wavefront 
stability better than 25nm result in speckles signifi- 
cantly weaker than the stellar light's residual. 

— Thanks to the rotation of the array, the "speckles" will 
be decorrelated from the planet's signal (this is not 
true in classical imaging with telescopes). This effect 
is described in §4.2. 

— Postprocessing of individual frames allows, to some ex- 
tend, to identify wavefront errors (§4.2). 

— We propose to use telescopes of modest size (3-meter 
diameter), which will reduce the number of actuators 
(and their stroke) needed to actively control the wave- 
fronts. This will keep the amplitude of wavefront vari- 
ations small. 

This photon noise is left after subtraction of the cen- 
tral star light on each snapshot and alters the estimation 
of the (u,v) components of the image and creates noise 
in the final reconstructed image. This noise appears as 
speckles with a maximum spatial frequency equal to the 
maximum spatial frequency of the signal. Figure 11 illus- 
trates this effect: this simulated image of the solar system 
at 10 pc is to be compared with Fig. 8. In this example, 
the simulator was used to generate a sequence of 100 noisy 
snapshots of the solar system and a sequence of 100 noisy 
snapshots of the sun. Between each snapshot, the array 
was rotated by of a full rotation, resulting in a rota- 
tion of half a turn between the first and the last snapshot 
of the sequence. Each snapshot is a 972.s/j,m exposure, so 
that the total integration time is 27/?,/im. Each of the two 
images (Solar system and Sun) was reconstructed using 
the technique presented in §2.4. The techniques that can 
be used to reconstruct an image of the star without the 
planets, using the same set of snapshots, are presented 
in section §4.2. The difference between the two images is 
shown on Fig. 11. 

The Earth is the brightest point in this image, and 
the S/N of the Earth detection, as measured on this re- 
constructed image, is 10. Jupiter is easily visible. Mercury 
and Venus (S/N = 3) are seen at the detection limit, but 
Mars is lost in the noise {S/N < 2). In this particular 
example, the exposure time needed to reach S/N = 10 
for the detection of the Earth in snapshot mode is, from 
equation B.23, 812 s./xm. The non-redundancy of the ar- 
ray is responsible for a multiplication by 4.7 of this expo- 
sure time, while the absorption of the Earth's light by the 
coronagraph in some of the snapshots is responsible for 
another multiplication by 3.67. In this simulation, photon 
noise was also simulated to generate the reference image of 
the Sun without planets. This multiplies the noise in the 
final subtracted image by v^, and therefore requires an 




Fig. 11. Simulated reconstructed image of the solar sys- 
tem at 10 pc. Photon noise included. Total exposure time 
{Tgff) is 27 h (with a l/xm bandwidth). 



exposure time twice as long to reach the same S/N. When 
all these factors are taken into account, the expected ex- 
posure time to reach S/N = 10 is 7M7mn for a l/xm 
bandwidth (100% efficiency). 

As can be seen in the image, most of the flux of the 
Earth is not in the narrow diffraction peak, but is diffused 
in a wide "cloud" : this effect is due to the coronagraph 
and the use of a reconstruction algorithm that does not 
take it into account (lack of deconvolution) . This explains 
why the exposure time required to reach S/N = 10 is 
3.5 times longer than predicted above. Because this effect 
is predictable, a proper reconstruction alogrithm (decon- 
volution) could recover part of this loss. In this exam- 
ple, the combined side-effects of the coronagraph (lower 
transmission of the planet's photons and the use of an im- 
age reconstruction algorithm which does not take it into 
account) are responsible for an increase of the exposure 
time by a factor 12.8. Although this factor is high, it is a 
small loss compared to the gain that coronagraphy brings: 
without the coronagraph, the exposure time would have 
to be 190days [Ijim bandwidth) to reach the same S/N. 
Moreover, some of this loss can be recovered by using a re- 
construction algorithm that takes into account the corona- 
graph, and this loss is decreasing rapidly with the number 
of apertures. 

5.4. Spectroscopy 

In this concept, the light of the companion is spread in 
many diffraction peaks on the focal plane detc;c;tor for each 
snapshot exposure. To obtain a spectrum of the compan- 
ion, we have identified 2 possible strategies: 
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Use of a wavelength-sensitive detector array 

If each pixel of the detector array used in the focal plane 
of our concept is sensitive to the wavelength of the in- 
coming photons, it is possible to process the series of 
frames independently in wavelength bins. Spectral infor- 
mation is then accessible for each source in the field. 
Arrays of Superconducting Tunnel Junctions (STJs) seem 
to be the ideal detector for this purpose. STJs detec- 
tors are photon-counting devices, and are used to mea- 
sure the wavelengths of individual photons. An experi- 
ments ( Peacock et al. 199^ ) has demonstrated their abil- 
ity to measure the wavelengths of photons from the UV to 
the NIR, with a typical spectral resolution of 10 (spectral 
resolution decreases with wavelength). A six by six pix- 
els array of S TJs has already be e n used for astronomica l 
observations ( Rando et al. 20do| |de Bruijne et al. 2002| ). 
Development of larger arrays at longer wavelengths with 
a spectral resolution of about 30 at lO^m would offer our 
concept the ideal focal detector for spectroscopic imaging. 
It is however presently unclear if this technology will be 
used at such long wavelength, because of the requirement 
for a very low energy gap supraconductor. 

Dispersing the companion's Ught on a detector 
array 

Spectroscopy of a point source in the field could be 
achieved by spatial selection of the source's light in the 
focal plane. If the position of the source in the sky is pre- 
viously known, the position of the diffraction peaks of its 
image in the focal plane can be computed. A mask can 
then select the source's photons and they would then be 
spectrally dispersed on a ID array. The mask would be a 
regular 2D set of holes if the array is redundant. 

A possibly easier way to achieve the same result would 
be to produce an image without redilution of the pupil to 
concentrate the source's photons in one single diffraction 
peak. The light in this diffraction peak would then be spec- 
trally dispersed on a ID detector. In either case, it would 
be more efficient to simultaneously acquire on-source and 
off-source spectra to correct for variation of the null and 
of the background. If the array' sub-apertures are along a 
single axis (in the entrance pupil and exit pupil), the light 
in the focal plane can be concentrated in a narrow "slit" 
(by using cylindrical optics) and then dispersed on a 2D 
detector array. 

Because the exposure time required to acquire a spec- 
trum is significantly longer that the exposure time re- 
quired to image the planet, we assume that a spectrum is 
being recorded for a previously imaged planet. Since the 
position of the planet is known, it is preferable to keep 
the array in a rotation angle chosen so that the planet 
is not significantly affected by the coronagraph. We can 
thereferore use the result of Appendix B to compute the 
sensitivity of the concept in spectroscopic mode. With a 
0.1 overall efficiency (including the loss of efficiency due to 
the phase mask, discussed in §5.2), a spectrum (S/N = 3) 
with a resolution of 30 can be obtained in 45 minutes for a 
redundant array and about 3.5 hours for a non-redundant 



array. With the redundant array, a spectrum with a reso- 
lution of 30 and a S/N of 30 can be obtained in 3 days. 

5.5. Discussion 

The S /N ratio in this concept is dominated by the effect of 
the residual star light Nsti and the zodiacal cloud compo- 
nent Nz in the example studied. As can be seen on fig. 12, 
the relative contributions of those two terms change with 
I, the distance to the system. When the system is close- by, 
the main source of noise is the residual star light which 
leaks through the coronagraph because of the large an- 
gular diameter of the stellar disk. This leakage decreases 
rapidly with I {Ngti oc while is independent of 
I. For I > 14pc, the zodiacal component dominates the 
noise. Nz is only a function of the local zodiacal light back- 
ground and the wavelength (Equation B.13), while A^^;, 
Nsti and Nst2 increase linearly with the total collecting 
surface of the interferometer (all other things kept equal). 
For I > lApc, the sensitivity of the concept is driven by the 
ratio Npi/Nz, which can only be increased by increasing 
the total collecting surface of the array (increasing N or 
d). For I < 14pc, decreasing the baseline B would improve 
the detection performance, but at the expense of reduc- 
ing the maximum I for which a system can be detected 
(about 30 pc with B — 60 m). It is interesting to note that 
for / > 25pc, Nst2 > Nsti'. the angular diameter of the star 
is becoming small enough for the wavefront errors (A/50 
at 0.5 /im in this example) to be the largest source of leak- 
age. If the total collecting area is increased, then this could 
become the largest source of noise for I > 25 pc. The ex- 
ozodiacal light component is much lower than either Nsti 
or Nz at all values of L In this example, it would signif- 
icantly lower the S/N only if the exozodiacal component 
is more than 20 Zodis. At such concentrations of exozodi- 
acal dust, the spatial structure (clumps, arcs etc..) of the 
cloud would be a serious issue. This important conclusion 
is valid regardless of the type of interferometer/null : this 
is a general law that applies to all concepts of total collect- 
ing surface 42 : Nz dominates N^z for cloud contents 
less than about 20 Zodis. 



5.6. Comparison with other proposed techniques 

Other proposed techniques use monopixel detectors com- 
bined with time-dependent optical pathlength differences. 
The flux received by the pixel is the product of the 
brightness map of the object (Star-l-planet(s)-l-zodiacal 
light -|-exozodiacal light) by a "transmission map" deter- 
mined by the geometry of the interferometer. The trans- 
mission map is modified (rotation for example) so that 
the light of the companion can be extracted. Chopping 
can be used to improve the sensitivity by increasing the 
frequency of the planet's signal on the pixel. This syn- 
chronous monopixel detection requires a stable and deep 
null. It is also sensitive to spatial anisotropy of the exozo- 
diacal cloud, because of the poorer imaging capability. In 



Olivier Guyon & Frangois Roddier: A Nulling wide-field imager for 



11 



Signals inside the companion's PSF 



- 




Npl+Nsl 1+N.st2+Nz+Nez 


- 








Nz 






... Nstl 
Nsl2 " ""^----...^ 






-- ^ Nez 















distance (pc) 

Fig. 12. Contribution {in ph. s^^) of the companion (Npl), 
the star leakage (Nstl, Nst2), the zodiacal light (Nz) and 
the exozodiacal light (Nez) to the photon count inside the 
companion's diffraction peaks, as a function of the dis- 
tance of the observed planetary system. 



our imaging concept, on each snapshot, the "on-source" 
and "off-source" channels are simultaneously measured 
(on different pixels of the detector). It is therefore eas- 
ier to correct for varying background level and varying 
null efficiency. For example, a variation in nulling perfor- 
mance from frame to frame can be computed and cor- 
rected for by analysis of the snapshots: its signature on 
the detector is different from the planet's PSF. In our 
concept, we combine nulling and spatial decoupling of the 
planet's photons and the star's photons on the detector, 
whereas other concepts only use one of those 2 techniques 
(nulling for interferometers and spatial decoupling for sin- 
gle aperture imagers). In this work, we have demonstrated 
how combining the 2 techniques lowers the technical re- 
quirements characteristic of each of the techniques: 

(1) The central null does not need to be in 9^ ot 9^. A 9^ 
null is sufficient. 

(2) The tolerance for wavefront errors and cophasing er- 
rors are lower: A/20 at 0.5 /im for detection of an Earth 
at 10 pc. 

In most nulling interferometer concepts, the null relies 
on a particular geometry of the array. In this concept, 
however, the null can be achieved with any array config- 
uration, thus offering the possibility of a very good (u,v) 
plane coverage and excellent wide field imaging capabili- 
ties. Complex exozodiacal structures can be imaged and 
the risk of confusion between a planet and such struc- 
tures is much lower. Imaging concepts using pupil densifi- 
cation and phase mask coronagraphy have been proposed 
by Boccaletti et al. (Boccaletti et al. 2000) and Guyon 
& Roddier (Guyon & Roddier 200C). In the concept pre- 




sented by Boccaletti et al., there is no redilution of the 
pupil, and the authors suggest the use of a high number 
of apertures {N = 36) to obtain a clean field of view large 
enough for exoplanet detection. In the Guyon & Roddier 
concept, a small number of apertures {N = 6) is used 



Fig. 13. Schematic of a possible optical layout for the con- 
cept. 

along with pupil redilution for wide field of view imaging. 
In both concepts, the imaging capabilities were very lim- 
ited, and aperture synthesis was not explored. The S/N 
computations in Appendix B also apply to those two con- 
cepts. 

5.7. Implementation - Complexity 

Unlike most nulling beam-combiners, the number of opti- 
cal elements in our Fizeau imager does not increase rapidly 
with the number of sub- apertures. Figure 13 shows a pos- 
sible optical layout for this concept: in this figure, there 
are a total of 6 x -I- 3 reflecting surfaces and one trans- 
missive element (the apodization mask, which could also 
be made reflective), excluding delay lines and at least one 
beam splitter (probably in the densifled pupil) for wave- 
front sensing and fringe tracking. Ml, the primary mirror, 
focuses the light on M2 which, with MS, is used to image 
the telescope pupil on the apodization mask. M3 also col- 
limates the beam. Optics of the delay lines have not been 
represented in this simple figure. The M4 mirrors pick up 
the beams to create the densified pupil. The M5 parabolic 
mirror forms an image on the phase mask and recoUimates 
the beam before pupil redilution. The M6 and M7 mirrors 
redilute the pupil before M8, the final imaging mirror. The 
nulling and beam combination optics are very simple and 
only have 2 x A^ 4- 2 refiecting surfaces and one transmis- 
sive mask. For A^ > 4 This number of optical elements is 
lower than for a non-coronagraphic nuUer, for which the 
number of beam splitters required is about A^^/2. This 
solution seems especially attractive for large numbers of 
apertures. 

6. Other capabilities 

6.1. Non-coronagraphic imaging at 10 /um 

When the phase mask coronagraph is removed from the 
focal plane between the 2 densified pupils, the interfer- 
ometer is an aperture synthesis imager. With six 2-meter 
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apertures in a non-redundant configuration, this concept 
can have a full (u,v) plane coverage (up to 60m baseline) 
with half a rotation of the array. The field of view of the 
reconstructed image is only limited by the optics and de- 
tector, and can be several times \/d. With lOmas pixels, a 
2048 X 2048 detector array would offer a 20 " field of view. 
The angular resolution at 10 iim would be 35 mas. No 
modification of the optics, other than moving the phase 
mask out of the focal plane of the densified pupil, would 
be required. 



6.2. Imaging in the visible or Near-IR 

The wavefront errors acceptable for imaging (A/10) are 
much larger than for nulling (w A/500). Achieving the 
wavefront correction required for nulling at 10 al- 
lows the interferometer to do imaging down to visible 
wavelengths. In this concept, most of the visible photons 
would probably be used for wavefront sensing and aper- 
ture cophasing, and the reflectivity of IR-optimized optics 
might not be very high in the visible. For those reasons, 
the sensitivity of the interferometer would probably be 
much below the theoretical sensitivity of a 42 w? collect- 
ing area (six 3 m diameter apertures), but the angular 
resolution at 0.5 /xm would be less than 2 mas (60 to base- 
line) . A dedicated visible detector array would be needed 
to take advantage of this potential. With a 10k x 10k array 
with 1 mas per pixel, the field of view would be 10". The 
scientific interest of non-coronagraphic imaging would be 
greater for objects significantly fainter, at visible wave- 
length, than close-by [d < 20 pc) stars. For example, the 
apparent visible luminosity of bright AGNs and quasars 
is typically 100 times less that a Sun-like star at 20 pc. 
This reduced number of photons for wavefront sensing and 
apertures cophasing could increase wavefront /cophasing 
errors and bring the shortest possible imaging wavelength 
form the visible to the Near-Infrared. 



6.3. Off-source imaging 

The very wide field of view theoretically achievable by 
this concept makes it possible to have both a wavefront 
sensing source and a science target in the same field of 
view. It is possible to obtain a reasonable sky coverage 
without relying on dual-beam intcrferometry, which would 
considerably increase the complexity of the interferometer. 
With a limiting visible magnitude of 15 (similar to the 
limiting magnitude of adaptive optics systems on 3-meter 
telescopes), and an off-axis guiding range of up to 20", the 
sky coverage is 10 %. If the science target and the guide 
star are separated by more than X/d (1 " at 10 /im), the 
"mixing" of the object photons and the guide star photons 
is very low. If it were an issue, the nulling coronagraph 
could be used to further reduce this effect. 



7. Conclusion 

The phase mask nulling stellar coronagraph can be ap- 
plied to interferometric arrays with a small number of 
telescopes to produce clean wide field of view images. 
Rotational aperture synthesis allows full (u,v) plane cov- 
erage up to a maximum baseline of 60 m with only six 
2m-apertures. The use of an imaging coronagraphic tech- 
nique is a very powerful tool to detect and image Earth- 
size planets around other stars. Exposure times of a few 
hours are sufficient to image an Earth-like planet around 
a Sun-like star. A planet is easily distinguished from an 
exozodiacal disk thanks to this imaging technique. The 
sensitivity to background emission (thermal emission of 
the optics, zodiacal light) and residual star light is very 
good because of the high spatial selectivity of the image. 
This high spatial selectivity also allows good spectroscopic 
sensitivity, which could be fully exploited with the use of 
ST J detector arrays. The wide field of view and good imag- 
ing capability of this concept also make it a very powerful 
instrument for general astrophysics. This solution is very 
attractive for simplicity of the optical layout especially for 
large number (more than 10) of apertures. 
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Appendix A: Performance of the phase mask 
coronagraph 

A.l. The theory of the phase mask coronagraph for a 
circular aperture 



The phase mask c oronagraph ( Roddier fc Roddier 1997 , 
[Guyon et al. 1999 ) uses a small mask to dephase the light 
in the central part of the Airy pattern by half a period. 
The mask is a disk covering 43% the diameter of the first 
dark Airy ring. As shown on Fig. A.l, in the focal plane 
of the telescope, the light complex amplitude is multiplied 
by -1 inside the phase mask (phase shift of half a period) 
while it remains unchanged outside the mask. 

With a and j3 the position vectors in the focal and 
pupil planes respectively, and A(/3), -B(a), C(/3) and 
D{q.) the light amplitude functions in planes A, B, C and 
D of Fig. A.l, 

Bioc) = J7iA{(3)), (A.l) 

where 3^7 denotes the Fourier transform (and 77~^ the 
inverse Fourier Transform). For a circular entrance pupil 
of radius R, A{f3) = 1 if |/3| < R, and \B{a)\'^ is the Airy 
function. When the phase mask is introduced : 

B{a) ^ Bo{a) - Bi{a). (A.2) 

Bo{a) and Bi(cx) are the complex amplitudes outside and 
inside the area covered by the phase mask, respectively. 



Olivier Guyon & Frangois Roddier: 



A Nulling wide-field imager for 



13 




entrance 
pupil 



imaging 
focal plane 



Focal Plane B 




B, 



Bq- b, 



No phase mask outside the 

phase mask 



inside the 
phase mask 



inside+outside 



C + C 




-0 

Pupil Plane C 



Fig. A.l. The phase mask coronagraph for a circular 
aperture. 

The effect of the phase mask is to multiply -61(0;) by -1. 
With the same notation: 



C{(3) = J7-\B{a)) = Co(/3) - Ci(/3), 



(A.3) 
(A.4) 
(A.5) 



Increasing the phase mask radius increases Ci (/3) and de- 
creases Co(/9) inside the pupil. For a critical value of the 
phase mask radius, C(/3) is then brought close to zero in- 
side the pupil, as shown on Fig. A.l. A "Lyot stop" that 
blocks the light outside the pupil in plane C is used to 
reject the light diffracted outside the pupil. For an off- 
axis point source, the phase mask in focal plane B has 
no effect (the Airy pattern is far from the mask) and all 
the light is inside the pupil in plane C: there is very little 
extinction for an off-axis source. With this simple phase 
mask coronagraph design, the maximum extinction fac- 
tor (monochromatic on-axis point source, optimal phase 
mask size) for a circular pupil is 160. This simple design 
has been successfully tested in an experiment we carried 
out on an optical bench with a monochromatic light source 
(Guyon et. al., 1999). 

A. 2. Entrance pupil apodization 

When the phase mask size is optimal, the extinction factor 
for an on-axis point source cannot be better than 160 for 
a uniformly lit circular pupil. This can be seen in Fig. A.l: 
C{f3) is not flat inside the pupil but has a positive cur- 
vature on this representation. If the entrance pupil, A{(3) 
has a negative curvature (transmission increases towards 
the center of the pupil), the positive curvature of C(/3) 




Fig. A. 2. Apodized entrance pupils for the phase mask 
coronagraph : single circular pupil (left) and interferome- 
ter's densified pupil (right). 



can be canceled. By a simultaneous optimization of the 
transmission map of the entrance pupil and the size of the 
phase mask, it is possible to reach a total extinction for 
an on-axis point source. 

We have developed an algorithm that computes the 
optimal transmission map for the pupil by changing A{(3) 
according to the C(/3) residual inside the pupil: three it- 
erations bring the theoretical extinction factor from 160 
to 10^ (monochromatic light, on-axis point source). For 
a circular entrance pupil, the transmission map has a 
49% minimal transmission at the edges of the pupil and 
a total integrated transmission of 73%. The apodization 
mask for such a pupil can be seen on Fig. A. 2 (left). The 
phase mask size has to be increased when the apodiza- 
tion mask is used. Such apodization masks have been 
studied in great detail in the case of rectangular pupils 



( Aime et al. 2001 , Aimc et al. 2002 ). In the following dis- 
cussion, we will consider an apodized entrance pupil. 

The algorithm used to compute the apodization mask 
for a single aperture telescope can also be used in a multi- 
aperture pupil. The result of this algorithm for our 6 tele- 
scope interferometer densified pupil is shown on Fig. A. 2 
(right). It is interesting to note that even for an unfilled 
pupil, this apodization mask yields a infinite theoretical 
extinction for an on-axis point source. However, as the 
filling factor of the entrance pupil goes down, the total in- 
tegrated transmission of the apodization mask decreases. 
For our 6 telescope interferometer densified pupil, the min- 
imum transmission is 23% on the outside edges of the den- 
sified pupil and the total integrated transmission is 57.3%. 
Had the sub-pupils been hexagonal, the total transmission 
would have been 66% thanks to a better filling factor. In 
this study, we will consider an ideal nulling coronagraph in 
which the chromatism problem has been solved, equipped 
with an apodization mask in the entrance pupil plane. 

In this appendix, we also estimate the amount of resid- 
ual light in the coronagraphic focal plane arising from dif- 
ferent errors. 
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A. 3. The phase mask thickness 



A. 5. The spectral bandwidth 



The relation between the phase shift p and the thickness 
e of the mask is given by 



p={n-l) 



e2n 



(A.6) 



n being the index of the material used to make the 
mask. A phase mask thickness error of Se yields a phase 
error 

6e 2tt 



5p = {n — 1)- 



A 



(A.7) 



For an on-axis point source, and an apodizcd entrance 
pupil, the residual light amplitude distribution in the pupil 
plane is the Fourier transform of the complex amplitude 
error in the focal plane. 

Res{l3) = 5T[(n - l)Se x Bi{cx)] (A.8) 
The corresponding light contribution in the focal plane is 



The phase mask needs to be achromatic : the phase shift 
should be half a period over the range of wavelength used 
for the imaging. Such an achromatic phase mask could 
be realized by stacking several layers of different materi- 
als, whose index-wavelength curves have to be carefully 
chosen. Let's consider an observation integrating wave- 
lengths from Ao — 6X to Aq -|- 6X (flat spectrum) with a 
phase mask coronagraph optimized for Aq. Without an 
achromatic phase mask, the phase shift introduced by the 
mask is a linear function of the wavelength, and, by inte- 
gration of equation A. 11. over the bandwidth, we obtain 



(A.16) 



r 5e f iengtn is equivalent to 

L = / Res{f3fdf3 = ir^-fx / {TT[Bi{cx)]fd^{J^S^ the 2 are linked by 

J\)3\<R ^ J\0\<R 



where F3 is the light "leakage" arising from the wavelength 
dependance of the phase shift. 

Another source of chromatism is the variation of the 
Airy pattern size with wavelength. The change of wave- 
length is equivalent to a change of the phase mask size 



'\P\<R 



The fraction of the total incoming light of the point source 
that "leaks" into the final focal plane is 



^A 
A^ 



6p 



(A.17) 



^ (-) 



3^7[Bi(a 



In the apodized pupil scheme, J7[Bo{a 
inside inside the pupil. 



(A.10) 





Prom integration of equation A. 14, the light "leakage" is 
then 



F4 



J {TI[Bo{cx) + Bi{cx)]yda 



(A.ll) 



Numerically, 



With (^) = 0.001, Fi = 2.5 lO"*^. 

A. 4. The phase mask diameter 

When the mask diameter is changed from po to po 
the same approach leads to 



Fa 



0.567 X 

Ao 



-X3(^)'(A.18) 



(A.19) 



Res{l3) = !FJ[2B{oc)p,^\c,\<po+Sp] 



dp, 
(A.12) 



And the fraction of the total incoming light that leaks into 
the the final focal plane image is 



F2 = 



U<R mm^)po<\c.\<po+Sp])^dcX 



J {TI[Bo{cx) + B^iaWdcx 
for small values of 6p 

^ pI X !\c.\<R (3T[2B(a)5(|a| - po)]fda 



(A.13) 



x(^)2(A.14) 



where 5{\a.\ 
Numerically 



/(5T[So(a)+Bi(a)])2da po 
— Po) is one if |a| = po and zero elsewhere. 



1.702 X (^)2 
Po 



(A.15) 



For a 1 p,m bandwidth centered on 10 p,m, 0.15 percent of 
the light of a point source is "leaking" into the final im- 
age. This problem can be solved at the expense of adding a 
wavelength-dependent magnification assembly. Such a de- 
vice, made out of a few chromatic lenses, has already been 
successfully used to record fringes in white light (Roddier 
et. al.; 1980) and for speckle imaging (Boccaletti et. al. 
1998). This is will extend the bandwidth by a factor of at 
least 100, bringing down the light leak to less than 10~^ 
for a 1 /im bandwidth at 10 /im. Another very interest- 
ing solution to solve the chomatism problem, as suggested 
by Antoine Labeyrie (private communication), is to use a 
Bragg hologram as a phase mask. The mask could then 
be made so that each wavelength sees it as being the right 
size and having the right phase shift. 

A.6. Pointing errors and wavefront errors 

For a pointing error of Sa, the light "leakage" is, for small 
values of Sa 



With (^) 



0.001, F2 = 1.710- 



^5 = 



I\c\<R (^^[/ - Po)cos{e - esac)d0])^da 



J {J7[Boia) + Bi{a)])^dcx 



x|(5a|2(A.20) 
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Numerical simulations give, for small values of Sa. 
F, = 0.364 X 

Po 



(A.21) 



For higher-order phase errors, the computer simulations 
give 



F5 = 7 X (1 - Strehl) 



(A.22) 



where Strehl is the Strehl ratio of the PSF. This formula 
was tested for Zernikes polynomials 4 to 100 and 7 is then 
between 0.6 and 1, depending upon the Zernike polyno- 
mial. 

Appendix B: Point source detection sensitivity in 
a snapshot exposure 

In this appendix, we compute the point source detection 
sensitivity without making any assumptions about the 
image reconstruction algorithm. We consider a snapshot 
frame in which the distribution of the residual starlight is 
exactly known. The photons of the point source (planet) 
are concentrated in a series of diffraction peaks. We esti- 
mate the signal (planet's photons) and noise (photon noise 
from the residual starlight, background, zodiacal light and 
exozodiacal light) in these diffraction peaks. This compu- 
tation does not take into account the problem of confusion 
with other structures (exozodiacal light concentrations, 
other planets), a problem for which this concept is well 
suited. 

The computations are first done for a redundant array 
and we then give quantitative measures for the decrease of 
S /N due to the non- redundancy of the array. We suppose 
that the planet's position on the sky is not on a null of the 
coronagraph. In a real observation of a peviously unknown 
planet, there is a probability for the planet to be occulted 
by the coronagraph. This effect is quantified in §5.2. 

We consider a space interferometer of N telescopes of 
diameter d spread along a baseline B (end to end), oper- 
ating at A = 10 fim (spectral bandwidth of 1 fj,m). The 
observation target is a star of angular diameter R^t of 
light flux Fst (in ph.s^^ .m^ .fim~^) as seen from Earth 
and a companion with a light flux Fpi at an angular dis- 
tance a from the star. At 10/im, the ratio between Fst 
and Fpi for a Sun-Earth system is approximately 5 x 10^. 
We give analytical expressions for each of the terms of 
the S/N expression. We also give numerical values for 
the observation of a Sun-Earth system at 10 pc {Fst = 
5xl0^ph.s~^ .m~^ .^m^^ ,Fpi = lph.s~^ .m~^ .fim~^ ,Rst — 
2.26 x 10~^rad), with a 1 Zodi cloud (face-on). Whenever 
numerical values are given, we consider a 6 aperture in- 
terferometer, 3 m diameter each with a baseline of 60 m 
(A^ = 6, B = 60m, d = 3m). 

The companion's photon count 
The size of a diffraction peak in the PSF is X/B. There 
are {B/{d x ^/N)Y such diffraction peaks in the PSF (re- 
dundant array). With t the average transmission of the 
apodization mask in the pupil (0.57 for a 6 aperture inter- 
ferometer), r(Fp/7rd^)/4 is the total number of photons 



from each aperture in the diffraction peak of an aper- 
ture. Since the area of this diffraction peak is (A/d)^, 
T{Fpi-Kd'^) / {'iX^) is the peak intensity in the image of the 
companion by one telescope. Because N telescopes' images 
are combined coherently, the peak intensity of the recom- 
bined image is N"^ times the peak intensity of the com- 
panion image in one telescope. The photon count inside 
the diffraction peaks of the companion's PSF is obtained 
by multiplying the photon count per diffraction peak by 
the number of diffraction peaks: 



Fpiird^ 
4A2 



/ A , 2 , B 



N, 



pi 



N X T- 



Fpind^ 



(B.l) 



(B.2) 



Because we have made the assumption that the array is 
redundant, we obtain the expected result that all of the 
companion's light "collected" by the interferometer is in 
the diffraction peaks of the PSF. The numerical value for 
our example is 



pi 



24.17 ph. s~'^.fim~ 



(B.3) 



The residual light from the star 

Because the angular diameter of the Solar-type stars at 5 
to 20pc is 0.5 to 2mas, in this computation, we do not need 
to worry about the star's angular diameter being larger 
than the distance between consecutive diffraction peaks in 
the PSF of the interferometer. The star is also smaller than 
the angular resolution of the interferometer we consider. 
Therefore, for a given direction, the residual light from a 
point source at an angular distance 9 from the optical axis 
is proportional to 9^ (Appendix A). When using an array 
that has no strong preferential direction (unlike the (u,v) 
plane coverage optimized arrays), the null is symmetric 
and the off-axis residual flux fraction equation, is, as ob- 
tained from our simulations ( Guyon fc Roddier 2000| ), 

Fn = (1.6 X -)2 X 02 



A 



for 



< 



X 

1.6xB ' 



Fr - 1 
for 



> 



1.6xB ■ 



In this work, we will consider Rst to be smaller than 
^ By integrating this equation over the stellar disk. 



1.6xS • 



the residual light from the on-axis star is 



B i?2 
stres^Fst X (1 .6 X X . 



(B.4) 



The residual light from the star is mostly concentrated 
in a small ring around the optical axis, and the detection 
of a companion is thus easier if the angular separation is 
important. However, for this simple estimate, we will con- 
sider that this residual light is spread evenly inside the 
PSF envelope (X/d). This is a very conservative assump- 
tion that will yield to pessimistic values for the detection 
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of companions that are not in the bright ring of residual 
light from the star. The residual star light is spread over 
^, as is the light of the companion. The diffraction peaks 
from the PSF of the companion occupy of this area (re- 
dundant array). Therefore, from (B.4), the photon count 
due to this residual is 



The numerical value for our example is 
Nsti = 4721 ph. s^^.fim^^. 



(B.5) 
(B.6) 

(B.7) 



The wavefront errors 

With the same conservative assumption as for the residual 
light from the star estimate, and using (A. 21), with 7=1, 



Nst2=T^F,til~ Strehl). 



(B. 



In this formula, Strehl is the Strehl of the PSF in the 
densified pupil scheme, and is to be related to densified 
pupil wavefront errors, which include cophasing errors be- 
tween the sub-apertures. For small wavefront errors, the 
Strehl and the variance of the wavefront phase errors 
are related through 



Strehl = e^" 



(B.9) 



By inserting equation B.9 into equation B.8, for small 
wavefront errors, 



Nst2 = T — F,t<7^ 



(B.IO) 



Nsti — Nst2 for a = 1.53 lO^^rad, which corresponds to 
a A/410 precision on the wavefront at 10/im, or A/20 at 
0.5/im. Reaching this accuracy on the wavefront is possible 
with smooth polishing of all optics and good cophasing of 
the sub-apertures. In this study, we consider a wavefront 
flat to A/50 at 0.5^m, for which Nst2 = 0.16 x Nsti^ and 
we adopt 



Nst2 = 755 ph. s ^.fim ^. 



(B.ll) 



The instrumental thermal emission and zodia- 
cal light 

With Bz the background contribution from zodiacal 
light {phs~^m~^sr~^), the corresponding photon count 
[phs^^) inside the companion's PSF diffraction peaks is 



TT 

4 



Nz ^ X Bz X A^. 



(B.12) 



(B.13) 



We adopt a value of 15MJysr ^ for the zodiacal light 
component at lO^m. At lO/im, 



IMJysr^^ = 1.5 x 10^'^ phfim~^sr~^m~^s~^ 



(B.14) 



hence, 

Bz ^ 2.25 X lO^^phfim-^sr-^m-^s''^, 
and 



Nz = 1.01 lO^ph.nm 



-1 s-\ 



(B.15) 



(B.16) 



The thermal emission of the optics of the interferometer 
(telescopes -I- recombination optics + detector) can be in- 
cluded in Nz, but by cooling the optics below 40 K, this 
contribution is smaller than the zodiacal light. 

Exozodiacal light component 
Estimating the amount of exozodiacal light photons in 
the diffraction peaks of the companion's PSF is harder 
than what we did for the zodiacal light component be- 
cause the exozodiacal light cannot be considered as uni- 
form across the field of view of the interferometer. For a 
given exozodiacal cloud model, the photon count will de- 
pend upon the relative position of the companion and the 
star (center of the exozodiacal cloud). For this example, 
we adopt the following model for the exozodiacal cloud 
([Reach et al. 19951): 



— The cloud is seen face-on. 

— The cloud vertical optical depth decreases as r~° '^^, r 
being the distance from the star in Astronomical Units. 

— The grain temperature decreases as 

— The inner edge of the cloud is at r = 0.02 (in the solar 
system, distance at which the dust temperature is 1300 
K, the dust sublimation temperature). 

We adopt a value of 25 MJy/sr at r = 1 for a 1 Zodi 
disk, and we consider a z Zodi disk. With this model, the 
surface brightness of the exozodiacal cloud is, in MJy/sr, 



EZ(r) 



1117.53 X zxr 



-0.37 



o4.8xr0-42-i 



or, in p/i.s ^./im ^.m ^.sr ^, 

1.676 10^5 X z X 
EZ[r) = — 



o4.8xr0''2_i 



(B.17) 



(B.18) 



for r > 0.02, and for r < 0.02. In the final image, the 
companion light is concentrated in a regular pattern (re- 
dundant array) of diffraction peaks, separated by ^^^^ 
one from the other in each axis. There are also multi- 
ple images of the exozodiacal cloud, each separated by 
the same distance on each axis. The shift between 

these two regular patterns is given by the position of the 
companion relative to the star. Because of the periodicity 
of these patterns (redundant array), the number of exozo- 
diacal photons in the companion's PSF is a 2-D periodic 
function of the relative position of the companion and the 
star. By summing the exozodiacal cloud surface brightness 
over the number of diffraction peaks {-j^^) of the com- 
panion's PSF and the number of images of the exozodiacal 
cloud, we can compute Nez, the number of exozodiacal 
photons on the companion's PSF diffraction peaks 



Ez = Nt— X {—Y X Bez, 

4 ID 



(B.19) 



Olivier Guyon & Prangois Roddier: A Nulling wide-field imager for 



17 




Fig. B.l. (a) True exozodiacal light distribution in a 1 
Zodi system at 10 pc (spatial coordinates are in AU). (b), 
(c) and (d) Exozodiacal light Equivalent background level 
when this system is observed respectively with a 6, 10 and 
20 apertures interferometer (spatial coordinates in units of 
6). In (b), (c) and (d), the effect of the nulling coronagraph 
has been simulated by multiplying the Exozodiacal light 
Equivalent background level by the transmission profile of 
the phase mask coronagraph. 



where Bez is the Exozodiacal light Equivalent back- 
ground level (in p/i./xm~^.s~^.m~^.sr~^) and 



Bez = J2 EZ{^/{x + k X 5)^ + {y + I X 5)^), (B.20) 

k,l 



where k and I are integers which are used to compute 
the discrete points of the exozodiacal cloud wich will con- 
tribute to Nez and 6 is the distance between successive 
diffraction peaks, 



S = 



XxVN 
B 



Bez is a J-periodic function of x and y and is shown 
(for < a; < (5 and < 1/ < (5) in Fig. B.l ioi N = 6, N = 
10 and N = 20. As a comparison, the true exozodiacal 
cloud surface brightness is also shown in the same figure. 

This figure illustrates the effect of having a sparse 
aperture: the relatively small angular spacing between 
diffraction peaks (or fringes in 1 dimension) results in a 
"mixing" of a large amount of exozodiacal light into the 
light of the companion. In this figure, Bez was computed 
for a 60m baseline interferometer at lO/um observing a 1 
Zodi cloud at 20 pc. Bez is independent of the individual 
sub-apertures diameter. With 6 sub-apertures, the median 
Exozodiacal light Equivalent background level across the 
field is 603 MJy.sr-i and 5 is 0.8 AU (355 MJy.sr"! and 
1.1 AU for 10 apertures, 173 MJy.sr''^ and 1.5 AU for 
20 apertures). By replacing Bez by this value in equation 
B.19, we obtain (for 6 apertures, 3 m diameter each, 60 m 
baseline): 



NEz = Q0.5ph.s ^.fxm ^. 



(B.22) 



Signal to noise ratio 
The signal to noise ratio for the detection of a companion 
is (noiseless detector) 



(B-21) ^ = ^ 



pi 



'pi 



-N 



stl 



■ N,t2 +Nz + Nez ' 



(B.23) 
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where T^ff is the effective exposure time (in s-iim). With 
the values computed above, we obtain a S/N of 3 in 
lOl.s./im effective exposure time, corresponding to a 17mn 
exposure time with a 0.1 efficiency (transmission of optics 
multiplied by the quantum efficiency of the detector) and 
a 1/im bandwidth. The residual star light and the zodiacal 
light are by far the most important sources of noise in this 
computation. 

Non-redundant arrays 
When using a non-redundant array, the PSF of the com- 
panion is less contrasted: its light is spread over a larger 
area. With a redundant array, most of the light (light in- 
side the diffraction peaks) occupies a fraction Fred of the 
PSF, with 

Fred = ^- (B.24) 

For a non-redundant array, it is difficult to estimate which 
diffraction peaks should be considered as part of the the 

area in which the companion's signal should be taken into 
account for the S/N computation. Rather than measuring 
the signal in discrete domains, as was done for the re- 
dundant array above, one should deconvolve the snapshot 
image to recover the companion. However, such a decon- 
volution is not easy to achieve because of the position- 
dependency of the PSF introduced by the use of a corona- 
graph, and is beyond the scope of this work. For this sim- 
ple estimate, we model the non-redundancy of the array as 
an increase of the relative area F of the diffraction peaks 
in the PSF and a decrease of the flux they contain. The 
number of star leakage, zodiacal and exozodiacal photons 
"mixed" with the companion's photon is then increased by 
relative to the redundant array configuration. With 
a non-redundant array such as the one shown in Fig. 1, 
most of the companion's photons arc contained in half of 
the total PSF area: in the PSF of the non-redundant 6- 
apertures array, 80% of the flux is contained in 50% of the 
PSF's area. Therefore, with F = 0.5 we expect the noise 
term of equation (28) is to be multiplied by 




(B.25) 



which is equal to in our example, and the signal 
(planet's photons) is multiplied hy E = 0.8. Note that the 
photon noise from the companion itself is not increased, 
but since it represents a very small fraction of the noise, we 
can ignore it in this simple estimate. This simple estimate 
shows that using a non-redundant array will increase the 
exposure time required for detection of the companion by 
approximately j^^, which in our example is 4.7. With a 
larger number of apertures, the difference of S/N becomes 
larger between a redundant and a non-redundant array, 
but it also becomes easier to obtain a good (u,v) plane 
coverage without using a very non-redundant array. 
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